A method to find unstable periodic orbits for the diamagnetic Kepler problem 
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A method to determine the admissibility of symbolic sequences and to find the unstable periodic 
orbits corresponding to allowed symbolic sequences for the diamagnetic Kepler problem is proposed 
by using the ordering of stable and unstable manifolds. By investigating the unstable periodic 
orbits up to length 6, a one to one correspondence between the unstable periodic orbits and their 
corresponding symbolic sequences is shown under the system symmetry decomposition. 
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• Unstable periodic orbits (UPOs) as the skeleton of chaotic systems provide a powerful tool to analyze the property 
of low-dimensional dynamics As a simple physical system, the diamagnetic Kepler problem (DKP) displays 
I— I important roles of UPOs in the classical and quantum chaos (2|j^. For the chaotic system, the UPOs are filled in 
^ the phase space and their total number roughly increases with the periodic length. Many UPOs can be found by 
. searching along three immediate symmetry lines in the configuration space jj] . In view of the scattering of a particle 
• ' in the four-disk system, symbolic dynamics has been established and a method to find UPOs by contracting a 
. ^ rectangle in terms of the symbolic dynamics has been given . By considering the stretching and wrapping in the 
1^ lifted space, symbolic dynamics has been established Q]. In this paper, we will present a method to find UPOs based 
on the ordering of stable and unstable manifolds and describe a one to one correspondence between the UPOs and 
their corresponding symbolic sequences for the DKP. The paper is organized as follows. In Sec. O, we describe the 



— I ' method to locate initial values of UPOs in terms of symbolic dynamics. In Sec. Ill, by investigating the UPOs up to 
^ . length 6, a one to one correspondence between the UPOs and their corresponding symbolic sequences is shown under 



the system symmetry decomposition. Finally, in Sec. |rV|, some conclusions are given. 

II. METHOD TO LOCATE THE INITIAL VALUES OF UPOS 

The classical dynamics of a hydrogen atom with zero angular momentum in a uniform magnetic field B along the 
z-axis is described by the "pseudo" Hamiltonian (using the semi-parabolic coordinates): 



p : ^ I + I - ,(^2 ^ ^2) ^ 1^2^2(^2 ^ ^2) ^ 2, (1) 

' where e = E^'^^^^ is the scaled energy, h has C^y and time-reversal symmetry. 
H ] In view of the Birkhoff canonical coordinates, a Poincare section has been chosen along a counterclockwise contour 
0- We choose the /i and u coordinate axes as a Poincare section and determine a Poincare map in terms of the 
entering direction of a orbit. In the first quadrant, the arc coordinate of the map is defined as s = —fj, along the ^ axis 
and s = V along the v axis. Using the available transformation jrppj ~^ Sy we get s G [—1, 1). The second coordinate 

of the map is v — , ^'^ ^ along the u axis and v — , ^ along the v axis. According to the C4 symmetry, the 

Poincare maps in the second, third and fourth quadrants can be obtained by rotating the map in the first quadrant 
with angles 7r/2, tt and 37r/2, respectively, around the center anti-clockwisely and adding 2,4 and 6 to s. Thus, the 
dynamics on the Poincare surface is represented by the map on the annulus s S [—1,7) and w € [—1,1]. The reduced 
domain (s e [0, 2), v € [—1, 1]) is constructed by considering the rotational symmetry. In the same way, the minimal 
domain (MD) (s € [0, 1), v e [—1, 1]) is taken by discarding the 7r-rotation around the origin. 

To illustrate the method to find UPOs, we first discuss the ordering of stable and unstable manifolds in the Poincare 
section. In Fig. 1, the stable and unstable manifolds in the MD display the local and global constricting, stretching 
and folding processes of the dynamics. When the tangencies associated with folding of one family of the manifolds 
occur, a partition line through them need to be introduced. Using the method given in Rcf. [^, we roughly determine 
the partition line "Co in Fig. 1. The demarcation line between two regions with the different rotation numbers in the 
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lifted space is marked by •-Bo- In terms of the natural ordering in the lifted space and the tangencies of manifolds, we 
obtain the region partition in the MD with the symbols {Lq, Rq and Ri) and the symbolic ordering of single letter 
for forward and backward sequences 

• Lq < •Rq < •Rl, Lq» < Rq» < Rf . (2) 

Except the ordering, the parity of a leading string in the sequences is defined by the oddness of the total number of 
Rq. Only odd leading string reverses the ordering. 

After the region partition and symbolic ordering are introduced, the families of the stable and unstable manifolds 
make a curvilinear coordinates system in the MD. All points on a stable (unstable) manifold converge to each other 
when iterated forward (backward). In terms of the language of symbolic dynamics, the points in each stable (unstable) 
manifold have the same forward (backward) sequence, i.e., the same ordering. A family of stable (unstable) manifolds 
is ordered according to their transverse intersections with an unstable (stable) manifold. The ordering is well-defined 
as long as there are no tangencies between the two manifolds. The occurrence of tangencies changes the ordering of 
folded stable or unstable manifolds, i.e., two families of sub-manifolds divided by the partition line have the reverse 
ordering. In Fig. 1, the ordering of stable (unstable) manifolds increases from the left-bottom (left-top) to the right- 
top (right-bottom) along each unstable (stable) manifold. For every point in the MD, there exists a stable or unstable 
manifold, which passes through the point and joins with the boundary line •Dq- 

The ordering of stable and unstable manifolds in the Poincare section can be displayed by using a and /? in the 
symbolic plane. Every forward sequence •SqSi • • • s„ • • • may correspond to a number a, represented in base 3, between 
and 1. The correspondence of symbol Sj to number G {0, 1, 2} is given as Lq ^ 0, Rq ^ 1, Ri ^ 2 if the leading 
string sqSi ■ ■ ■ Sj_i is even, and as {Lq, i?o, Ri} — > {2, 1, 0} for the odd case. Here, the number a is defined as 



a = £63-(^+i). (3) 



i=0 

Similarly, the number f3 for a backward sequence • • • S-m ■ ■ ■ S-2 is defined by 

oo 

P = J2^jS~'' (4) 

i=i 

where rjj € {0, 1, 2} is determined by S-j and the oddness of the leading string S-^+i • • • s_2S-i- In this way, forward 
(backward) sequences are ordered according to their a (/?)-values. In particular, at the boundary line s = 0+, a (/?) 
increases (decreases) monotonically from f = — ltow = l. 

In Fig. 2, the symbolic plane has a one to one correspondence with the MD. Many point sets are separated by the 
boundary line •Dq, the partition line •Co and their forward or backward maps. For a given symbolic sequence, its 
ctexact and (3exact valucs may be calculated by using (3) and (4). In the moving process of a starting point, ttapprox 
and (iapprox of the current point are respectively determined by using its forward and backward sequences formed in 
50 forward and backward maps of the current point in the MD. Our method is described as follows: 

(i) We run near a^xact from v=-l in the boundary line 5 = 0"*" and then iterate v to aexact by using the half-division 
method. In the iteration process, if both aapprox and (iapprox for one of the two points are close to aexact and Pexact 
with a given error (10^^), we obtain an initial value of the UPO corresponding to the given symbolic sequence. The 
UPO is that passing through the origin. 

(ii) If \aapprox — OLexact\ < 10~^, i.e., the difference between the values aapprox of two points surrounding aexact 
is less than the error, we may make a stable manifold from the middle one between the two points. Along the 
stable manifold, aapprox remains approximately unchanged and Papprox increases monotonically. When (^approx is 
near (3exact, i-e., the current point leaps from one side of (3exact to another, one should decrease the step size to achieve 
\liapprox — Pexactl < 10~^. Howevcr, in this case, the requirement \aapprox — aexact] < is not met. We may 
make an unstable manifold from the middle point and choose a direction according to aapprox- Along the unstable 
manifold, Papprox remains approximately unchanged and aapprox is monotonically changed. At this time, if the 
change of Papprox or aapprox decreases monotonically, the above process will be continued until both the requirements 
\aapprox - aexact] < 10~^^ and \/3approx " Pexact\ < 10~^^ are met. Thus, in terms of the differences \aapprox - aexact] 
and \f3approx — Pexaci\, wc may extract an initial value of the UPO corresponding to the symbolic sequence, which 
is allowed. The UPO is that non-passing through the origin. If the change of Papprox or aapprox does not decrease 
monotonically in the 10 circulations, i.e., the algorithm does not converge, the above process will be stopped. Thus, 
we give up the process for aapprox and repeat (ii) for Papprox- 

(iii) If \aapprox — aexact \ > 10~* cvcn when the step of the half-division method is enough small, i.e., aapprox jumps 
between two nearly invariant values around aexact, we give up the process (ii) for aapprox and repeat it for papprox- 
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(iv) If IPapprox — Pexact\ > 10~^ cvcn whcn the step of the half-division method is enough small, i.e., Papprox also 
jumps between two nearly invariant values around Pexacti we may say the symbolic scqiience is possibly forbidden. 

Using the above method for the given symbolic sequence, we can obtain an initial point of its corresponding UPO 
or a possible estimation of forbiddenness. For the latter case, we may further repeat the above steps to determine 
the admissibility of symbolic sequences produced by shifting the given sequence in its length. We stop the process 
as soon as the symbolic sequence is allowed. If all symbolic sequences generated by shifting the given sequence in its 
length are possibly forbidden, we may say the given symbolic sequence is forbidden. Finally, we can determine the 
allowance or forbiddenness of the given symbolic sequence and extract an initial point of the UPO for the allowed 
symbolic sequence. 

To illustrate the main running process of the above method, we take the symbolic sequences B^Ri and LqRqRx 
as two examples. For the sequence (-R^i?i)°° • {RlRi)°° , using (3) and (4), we get aea^ac*=0.5384615384615387 and 

/3ea;act=0. 8461538461538460. The entire process to find the UPO is displayed in Fig. 1 and Fig. 2. First, we begin 
from the point (so7^'o)=(10~^°,-0. 99999) in Fig. 1 and take its symbolic sequence to calculate {a,fi) in Fig. 2. We run 
along the line s — so with a step Awq and determine the last point near (Xexact by comparing Uapprox with a^xact- 
Then, we use the half-division method to decrease the difference between Uapprox and aexact- When the difference 
is less than 10~^, we draw the stable manifold from the point in Fig. 1. Along the stable manifold, we compare 
Papprox with Pexact- When the difference between them is less than 10~", WC C0nipa,r6 Otapprox witll OLexact 

and get 

their difference 0.0006378404570228. The last point is marked by a in Fig. 3(a). According to the comparison, we 
draw the unstable manifold from the point. When \aapprox — ctexactl < 10~^^, we compare Papprox with pexact and 
get their difference 0.0000029487355183. The last point is marked by b in Fig. 3(b). Furthermore, we continue the 
above steps with the given error 10~^^. In the second circulation, when \l3approx — Pexactl < 10~^^, the difference 
\<^approx — ctexactl is 0.0000000000441224. The last point is marked by c in Fig. 3(c); when \aapprox — ctexact] < 10~^^, 
the difference \Papprox - (3exact\ is 0.0000000000000114. The last point is marked by d in Fig. 3(d). In the third 
circulation, when \l3approx — Pexactl < 10^^^, the difference \aapprox — oiexact\ is also less than 10~^^. Thus, we get 
the point (s, w)=(0.6872027097581007,-0.3670108804483804) marked by e in Fig. 3(d) corresponding to the sequence 
{RqRi)°° • (RqRi)^. Also, we may obtain its forward and backward sequences with 50 letters to conform the 
result. Finally, using Newton method, we exactly extract the UPO. For the sequence {LoRoRi)°° • {LoRoRi)°° , we 
get ae^act=0.1428571428571428 and /?ea:act=0.8571428571428572. First, beginning from the point (so,t;o)=(10"^° ,- 
0.99999), we run along the line s = sq with a step Avq and determine the last point near Uexact by comparing 
ctapprox with ttexact- Then, using the half-division method, we decrease the difference between Uapprox and aexact- 
When Avo < lO'^^, the Uapprox values of the last and current points are 0.1428571417083783 and 0.1428571738737871, 
respectively. Since Auapprox > 10~^, we give up the process for Uapprox and repeat it for Papprox- When Avq < 10~^^, 
the Papprox values of the last and current points are 0.8553113553113367 and 0.9046940713607194, respectively. Since 
Aaapprox > 10~*, the process for Papprox is stopped. Thus, the sequence LqRqRi is possibly forbidden. Shifting 
the sequence LqRoRi, we get RqRiLq and extract aea;act=0.4285714285714285 and /3ej,act=0.2857142857142856. We 
repeat the above process. When Auo < 10" ^■'^ for aapprox, Aaapprox = 0.0000000964962270 > 10"®. When Avq < 
10-15 for Papprox, APapprox = 0.0164609053497942 > 10"*. The sequence RqRiLo is possibly forbidden. Shifting 
the sequence RqRiLq, we get RiLqRo and extract aexact=0.7142857142857142 and /3exact=0.5714285714285717. We 
repeat the above process. When Ava < 10^^^ for aapprox, Aaapprox = 0.1481481481481481 > 10"*. When Avq < 
10-15 for Papprox, APapprox = 0.0054869684499316 > 10"^ The sequence RiLqRq is possibly forbidden. Since the 
three sequences generated by shifting LqRqRi are possibly forbidden, the sequence LqRoRi is forbidden. 



III. CORRESPONDENCE BETWEEN UPOS AND SYMBOLIC SEQUENCES 

Before using the method to find UPOs, we discuss qualitatively the bottom-left region forbidden by the boundary 
line "Do in Fig. 2. Let us focus on the orbits near the origin in the configuration space in Fig. 4. Since the orbits are 

nearly linear, we can obtain the relation between two continuous Poincarc maps. A current point near the origin is 
described as (sq, vq), its forward and backward mapping points can be written as {si,vi) and {s-i,V-i), respectively. 

(i) If Vq € (—1,0), an orbit between the points {so,Vo) and {si,vi) is nearly linear and has not turning points. We 
have 



S, = « -SOV^/VO > 0, 

vo {so-l)+ soy/l-v^ * ^' (5) 



Vl 



where vi € (0, 1). We can obtain the symbolic description Lq for the point (so,fo) and Rq or i?i for the point (si,wi) 
in Fig. 1. The symbolic string (Lq) is forbidden due to the geometry of the orbits. When the coordinate vq changes 
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from -1 to 0, the symbolic sequence for the point (sq, ^^o) changes from • LqRqR^ to Rf'Ro • LqR^ and a {(3) 
increases (decreases) monotonically from ^ (1) to i (i). Thus, we get the left half part of the boundary line •I?o in 
Fig.2. 

(ii) If uo G (0, 1), an orbit between the points (s_i,w_i) and (sqj^^o) is nearly linear and has not turning points. We 
have 



V-1 = 



"0' 



where w_i G (—1, 0). We can obtain the symbolic description Rq or Ri for the point (sq, vq) and Lq for the point (s_i, 
in Fig. 1. The symboHc string (Lq) is forbidden due to the geometry of orbits. When the coordinate wq changes 
from to 1, the symbolic sequence changes from R^Lq • RqR^ to Rf^RoLo • Rf^ and a (0) increases (decreases) 
monotonically from -i (i) to 1 (i). Thus, we get the right half part of the boimdary line •Z?o in Fig- 2. 

Therefore, combining (i) and (ii), the bottom- left region forbidden by the boundary line 'Dq is obtained due to the 
linearity of the orbits between two continuous maps near the origin. 

In case of three symbols with the above exclusion rule, the associated connectivity matrix is 
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Ro 
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Lo 
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Ro 
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Ri 
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1 



The total number of symbolic sequences of the length n is Nn = trT" , which can be decomposed into the number AI^ 
of primary symbolic sequences of length d dividing n: Nn — X]rf|n ^^^d- By the Mobius inversion we get 

n ^ a 

d\n 

where /i(l) = 1, /^(fi) = if n contains the square of a prime and ^{n) = (— 1)*^ if n contains k prime factors. All the 
calculated results are shown in Table I. 

In Ref. 10 we have shown that the symbolic dynamics in the MD removes the rotational and reflectional symmetries 
and only preserves the time-reversal symmetry T, which just reverses the original sequence. We make a sequence set 
with different lengths by combining the three symbols. From the set, the symbolic sequences generated by shifting 
each primary sequence are removed. If a symbolic sequence generated by the time-reversal symmetry is only shifts of 
the primary symbolic sequence, it is also abandoned. Finally, all symbolic sequences including Lq should be removed 
from the set. Thus, the number of the sequence set is equal to Mn obtained from the connectivity matrix, which is 
shown in Table I. Under the system symmetry decomposition, only the time-reversal symmetry is preserved in the 
sequence set, which can be used in cycle expansion |p^ . 

Using the method mentioned above, we determine the admissibility of a sequence set up to length 8 and extract 
approximate initial points of the UPOs corresponding to the allowed symbolic sequences for e = 0. The number of 
allowed symbolic sequences is presented in Table I. In the primary sequences, their admissibility as physical UPOs 
can be further determined in terms of pruning fronts corresponding to the tangencies between stable and unstable 
manifolds and bound line 0. For example, in Fig. 5, the points (a.^aij^, PmR'^^R^) and {oimR„R-^^R„, P%RoRiRo) exist 
in the allowed •i?Q region, and the point (a,^jfl2, P,r-^r2) exists in the allowed uRi region. So, the sequence 
(i?o^i)°° • {RlRi)°° is allowed. The point {a,LoRoRn P*LoRoRi) "iocs not exist in the allowed "Lq or •i?Q or uRi 
region, but in the forbidden region of boundary line •Do- So, the sequence {LoRoRi)°° • {LoRoRi)°° is forbidden. 
We use the method to determine the admissibility of primary sequences up to length 8. The result is consistent with 
the above one. 

For a real number, our Sun- workstation has significant digits with the 16's space. Due to the restriction of precision, 
the exact UPO up to length 6 can be further exacted by using the Newton method. In order to check a one to one 
correspondence between the UPOs and their corresponding symbolic sequences, we display the allowed symbolic 
sequences up to period 6 in Table II and their some corresponding UPOs in Fig. 6. In terms of the rotational 
and reflectional symmetries, the Orbit Period is the same or twice or quadruple the Sequence Period, which can 
be determined in the configuration space. In the 38 UPOs of Table II, except the pairs of UPOs (5) and (6), (11) 
and (15), (14) and (16), (25) and (29), (26) and (32), (27) and (33) shown in Fig. 6, other UPOs have the different 
configurations and cannot be produced from one to another in terms of the time reversal. For each one of the 26 UPOs, 
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the timc-rcvcrsal orbit is the same as that produced by the rotational symmetry, or reflectional symmetry, or their 
combination. So. the time-reversal symmetry T is degenerate, i.e., the UPO and its time reversal orbit correspond to 
the same symbolic sequences in Table II. For the pairs of UPOs (11) and (15), (14) and (16), (26) and (32), (27) and 
(33), they have also the different configurations, but the two UPOs of each pair have the same one. For each one of 
the two orbits, its time-reversal orbit is another one, which cannot be produced from the given UPO by the rotational 
symmetry, or reflectional symmetry, or their combination. So, the time-reversal symmetry T is non-degenerate. The 
two orbits with the non-degenerate T symmetry correspond to the different symbolic sequences in Table II. For the 
pairs of UPOs (5) and (6), (25) and (29) passing through the origin, they have also the different configurations, but 
the two UPOs of each pair have the same one even when the time reversal symmetry is considered. When the right 
limit s = 0+ of the UPOs is taken into account in the configuration space, the intersection points of the UPOs (5) and 
(29), (6) and (25) should be removed from the origin into the third and first quadrants, respectively. At this time, for 
the pair of UPOs (5) and (6) with the right limit, the two lines of the UPO in the center should be removed to pass 
through the axes > and v respectively; for the pairing UPOs (25) and (29) with the right limit, the two lines 
of the UPO in the center should be removed to pass through the axis v > {). The pairs of UPOs (5) and (6), (25) 
and (29) with the right limit have the similar configurations to the above UPOs no-passing through the origin and 
the same results can be concluded. Thus, a one to one correspondence between the UPOs and their corresponding 
symbolic sequences is shown under the system symmetry decomposition. 



IV. CONCLUSION 

According to the ordering of stable and unstable manifolds, wc propose a method to determine the admissibility of 
symbolic sequences and to find UPOs corresponding to allowed symbolic sequences for the diamagnetic Kepler problem. 
By searching the UPOs up to length 6, a one to one correspondence between the UPOs and their corresponding 
symbolic sequences has been shown under the system symmetry decomposition. 
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Table I. Numbers of the total sequences (TS), the primary sequences (PS) and the allowed sequences (AS). 
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TS 


Nn 


2 
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20 


56 


152 


416 


1136 


3104 


PS 


M„ 


2 
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6 


12 


30 


65 


162 


381 


AS 


e = 
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12 


17 


39 


78 
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Table II. The allowed sequences encoding UPOs for e = 0. By Oi we mean the i-th orbit according to the numbering in the first column. 
The Sequence Period and Orbit Period are defined by the length of symbolic sequences and periodic points of UPOs in the Poincare 
section, respectively. 



No. Sequence Period Sequence Orbit Period Non-degeneracy Symmetry 

~ i 4 

2 2 LoRi 4 

3 2 RoRi 4 

4 3 R^Ri 6 

5 4 LoRlRi 8 

6 4 LoRiRq 8 T0O5 

7 4 LoRiRoRi 4 

8 4 R^Ri 8 

9 4 flgJjf 8 

10 5 LoRoLoRl 20 

11 5 LaRif^Rl 20 

12 5 LoRiLoRi 10 

13 5 Lo-Ri-Ro-Ri 20 

14 5 LoRiRoRl 10 

15 5 LoRlRl 20 To On 

16 5 LoRlRoRi 10 T o O14 

17 5 Rt,Ri 10 

18 5 flgiif 20 

19 5 RlRiRoRi 20 

20 5 RlRl 10 

21 5 RoRiRoRl 10 

22 6 LqRqLqR\RqR\ 12 

23 6 LoRoLoRi 12 

24 6 LoRlLoRl 6 

25 6 LoRtRi 12 

26 6 Lo-Ro-Ri 6 

27 6 LoRlRl 12 

28 6 LoRiLoRiRoRi 12 

29 6 LoRiRo 12 T o O25 

30 6 LoRiRoRi 12 

31 6 LoRiRoRiRoRi 12 

32 6 Lo-Rfflo 6 T0O26 

33 6 LoRfRo 12 T o O27 

34 6 RqRi 12 

35 6 RoRi 12 

36 6 i?gJ?ii?oi?i 12 

37 6 R^Rl 12 

38 6 R^Rf 12 
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FIGURE CAPTION 

Fig. 1 Stable and unstable manifolds with two partition lines in the MD. Circles show the entire process to find the 
point corresponding to the sequence {B^Ri)°° • {RqRi)°° in the MD. 

Fig. 2 A symbolic plane corresponding to the MD. Circles show the entire process to find the point corresponding 
to the sequence {RqRi)°° • {RqRi)'^ in the symbolic plane. 

Fig. 3 (a) A local sketch in Fig. 1; (b) A local sketch in (a); (c) A local sketch in (b); (d) A local sketch in (c). 

Fig. 4 A chaotic orbit displaying its linearity near the origin. 

Fig. 5 Determining the allowance of sequence {RqRi)°° • {RqRi)°° and the forbiddenness of sequence (Lo-Ro-Ri)°° • 
{LoRoRi)°° in the symbolic plane by pruning fronts. 

Fig. 6 Some UPOs corresponding to the symbolic sequences in Table II. 
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